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We briefly review some common diffusion-limited reactions with emphasis on results 
for two-species reactions with anisotropic hopping. Our review also covers single-species 
reactions. The scope is that of providing reference and general discussion rather than de- 
tails of methods and results. Recent exact results for a two-species model with anisotropic 
hopping and with 'sticky' interaction of like particles, obtained by a novel method which 
allows exact solution of certain single-species and two-species reactions, are discussed. 
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1. Introduction 

There has been much interest in the studv. of nonequilibrium systems through repre- 
sentation by simple stochastic modelsJM In particular, low-dimensional diffusion- 
limited reactions are tractable and in many cases, exactly solvable, examples of 
nonequilibrium systems. Since the early works of Ovchinnikov and Zeldovich,H Tou- 
ssaint and Wilczekja and Torney and McConnellEl who pointed out the importance of 
fluctuations of particle densitvjn the kinetics of low-dimensional reactions, a host of 
results have been published B~x?3 Despite their simplicity, common diffusion-limited 
reactions can describe complex phenomena. However, only recentl y atte ntion has 
been turned to the effect of biased diffusion on these reactions.Q"Ej'E2l Here, we 
review some common diffusion-limited reactions with emphasis on the effects of bi- 
ased diffusion in two-species reactions. Single-species reactions will also be briefly 
reviewed. 

The asymptotic diffusion-limited reaction behavior is associated with diffusion 
being the slowest, 'rate-determining' process in the system. For dimensions lower 
than the upper critical dimension, a mean-field theory using the classical, Smolu- 
chowski-type rate equations — see Kang and Rednerllil — breaks down since it does 
not take into account effects of spatial inhomogeneities of particle concentration. 
Specifically, the classical mean-field theory predicts, in one dimension (Id), for 
the kinetics of single-species reactions — A + A — > A and A + A — > (inert) — 
an asymptotic decay of the concentration proportional to i _1 instead of the exact 
result i _1 / 2 .ErE3'E3"E3 That diffusion is the mechanism leading to anomalously slow 
kinetics can be illustrated for the annihilation reaction, A + A — > 0, on a one- 
dimensional (Id) lattice, as follows. After a time t, a particle will diffuse over a 
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distance of order I = VDt, where D is the particle diffusion constant. Let us 
assume that the concentration at time t is C(t), so that the mean distance between 
particles is n(t) — 1/C(t). In the limit n ^S> I it is reasonable to assume that particles 
diffuse freely, i.e., perform independent random walks, and there is no annihilation. 
In the opposite limit, n <C I, it is very unlikely for a particle to actually diffuse over 
a distance (, without being annihilated. Thus the variation of the particle density 
with time will be determined by the 'balance' of diffusion and annihilation achieved 
when n(t) » £, and we get C(t) ~ t' 1 ' 2 . 

Once two particles annihilate, there develops a local fluctuation in the density 
and correlation functions. The spatial extent of the disturbance is n(t) and therefore 
the time scale required to erase this local fluctuation by diffusion is n 2 /D. However, 
if n w £, the next reaction event locally will occur on the same time scale! Therefore, 
diffusion is ineffective in 'mixing' the particles and as a result the Id dynamics is 
truly fluctuation-dominated, non-mean-field. Similar heuristic arguments have been 
advanccdH for the c?-dimensional case. 

Recent research has been focused on dimensions lower than the critical dimen- 
sion, where the mean-field description breaks down. In particular, Id problems 
are usually more amenable Ixl mathem at i cal analysis; in several instances exact 



solutions have been reportedJlj'0'll3'lla"E3'ca'ca Apart from their intrinsic interest, 
exact results, when available, serve to test theoretical and computational results 
and methods which are also applied to more realistic problems. 

Low-dimensional reaction-diffusion models and related simple dynamical sys- 
tems actually describe experimental phenomena as diverse as electron-hole recom- 
bination in semiconductors, soliton-antisoliton dynamics in quasi-l<i systems, evolu- 
tion of co m peti ng species in biology, aerosol dynamics, star formation, and polymc- 
rization,Ej'E3o"Eil etc. Despite their inherent complexity and diversity such systems 
share common features — some reviewed below. These 'universal' properties can be 
studied in simple stochastic lattice-dynamics models. 

Experimental investigations have been performed on quasi-ld systems. 0BeEe! 
In such systems the hopping rate in the favorable direction is about I0 4 times higher 
than in the plane perpendicular to that direction. Therefore, the system cam he 
regarded as Id to a good approximation. Notably, hexagonal crystals of TMMCcJcj 
represent the best quasi-ld system available experimentally. Excitons propagate 
preferably along the Mn 2+ chains, as compared to motion across chains; the latter 
distance is nearly three times larger than the in-chain separation. The resulting 
hopping rate between adjacent chains is eight orders of magnitude smaller than the 
in-chain hopping rate. Typical values for the hopping rates and lifetime are 10 11 - 
10 12 s _1 and 740 /is, respectively. This large lifetime excludes natural decay of the 
quasiparticles, which would yield an exponential relaxation, so that the dynamics 
is dominated by the quasiparticle reactions. It is also a good approximation to 
consider the reaction time as instantaneous. Indeed, a pair of excitons reacts in 
about 100 fs, i.e., at least one order of magnitude smaller than the typical hopping 
time. Theoretical studies of the effects of finite probability of reactions have been 
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reported in the literatureSE! El Next, we focus on analytical results for single- and 
two-species reactions. 

2. Some common diffusion-limited reactions 

The kinetics of diffusing-particle single-species Id systems reacting via the processes 
of coalescence (A + A — > A) or annihilation (A + A — > 0) on particle encounters, 
is now well understood. For these two basic processes the concentration decays 
for large times according to C(t) ~ t~ x / 2 (in Id). In general, C(t) ~ t~ d / 2 for 
d < 2, with aloffaxit^T^^orrection factor at d = d c = 2, and C[t) ~ for 
d > 2J'll3'H-0>Hfia'El'El>Ea Here d c denotes the upper critical dimension. The 
steady state, assuming additional processe s, e.g . , par ticle creation or back-reaction, 



has also been studied by several authors .I13'E3'IH3'E2I'E3 In particular, RaczcJ obtained 
exact results by mapping interfaces of the kinetic Glauber-Ising model to particles 
in the A + A <-» system. Amar and Familye3 pointed out that the Ising-model to 
reaction-dynamics mapping should be used with caution owing to the presence of 
subtle correlations in the initial conditions for non-zero values of the magnetization. 
Multiparticle kA — > reactions^nd other variants of single-species reactions have 
also attracted attention.E3E3'E3'EJ'LJ For k > 3, a mean-field approximation applies 



in Id. For the borderline case of k = 3, Monte-Carlo simulations^ and refined 
mean-field treatmentLj have recently succeeded in demonstrating the logarithmic 
correction in the concentration decay. 

The kinetic behavior of two-species reactions is far richer than that of the 
single-species reactions. Three two-species models belonging to different universal- 
ity classes are discussed below. Although there exist other variants of two-species 
models,E3'E3'E3"EJ we restrict the present survey to these few examples. In the 
'standard' AB model particles hop isotropically on a d-dimensional lattice to any of 
their neighboring sites. On encounters, AB particle pairs annihilate while for AA 
and BB pairs the two particles 'bounce off' each other. This models a hard-core 
interaction whereby each site can be occupied by at most one particle. 

The 'standard' AB model shows the non-mean-field concentration decay C(t) ~ 
t~ d / 4 for d < 4 and the mean-field decav C(tl~ t^ 1 for d > d c = 4, for equal con- 



centrations of the two species.! 



-field decav C(t) 



For unequal (initial) concentrations, 



the concentration of the minority species decays exponentially as C(t) ~ exp f— Xy t) 



for d = 1, C{t) ~ exp (-At/lnt) for d = 2, and C{t) - exp (-Xt) for d > 3.1 
The seminal works by Bramson an d Leb owitz summarize the exactly known results 
for the d > 1 'standard' AB modelJsN 

A 'driven' Id version of the AB model was recently introduced by JanowskyB: 
particles are constrained to hop in one direction only. When combined with hard- 
core interactions, the biased-hopping dynamics yields a new universality class. Phe- 
nomenological considerationsQ'cl relate this reaction model to the asymmetric simple 
exclusion process, for which the noisy Burgers' equation represents the continuum 
limit. Numerical simulations and the above connection to the noisy Burgers' equa- 
tion suggest the concentration-decay exponent i, compared to i of the 'standard' 
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AB model, in Id. Regarding the spatial domain structure, there is only a qualita- 
tive agreement in the literatureJM including phenomenological considerations and 
numerical simulations. 



3. Exactly solvable 'sticky' two-species annihilation model 

A third_.]zaxiant of the AB model in Id has been introduced recently by Kra- 
pivsky.E3'Eil'E3 In this system same-species particles interact by sticking together. 
Same-species particle clusters coalesce on encounters — no hard-core interaction is 
present — and they form a larger cluster with size given by the total number of par- 
ticles in the merging clusters. When two clusters of different-species particles meet 
AB-annihilation takes place. The resulting cluster has particle number given by 
the particle-number difference of the reacting clusters, while the species is that of 
the 'parent' cluster with the larger number of particles. Numerical simulations and 
scalingEJ suggest d c = 2. 

Synchronous-updating, cellular-automaton-like, 'sticky-particle' lattice dynam- 
ics models have been considered recently. Exact results have been obtained for 
'sticky' diffusion-limited single-speciesj£3 and single- and two-speciesHEil reactions. 
They have been treated by a novel approach which yields a unified exact solution of 
the coagulation and annihilation single-species reactions and the 'sticky' two-species 
reaction defined above J13 The approach comhineaand extends several techniques de- 
velopedin_charge-coagulation model studieaZSrO and earlier single-species reaction 
studiesEil'E3 of cellular-automaton-type models. A combinatorial argument is then 
utilized to relate thejcluster history' at time t to the set of compatible initial config- 
urations (at t — 0)H Exact solvability for the two-species reaction is based on the 
linearity of certain dynamical evolution equations, supplemented bv a mathematical 
development yielding a closed form for a non-trivial double-sum£3 

For the single-species reactions the expected power-law decay with exponent 
^ was confirmed, regardless of the hopping rate bias. Moreover, for large times 
(only) the effect of the hopping-rate anisotropy can be fully absorbed in the time 
dependence, i.e., it only changes the value of the diffusion constant. The method also 
allows one to extend the exact relation, 2p) = 2Cq,(t;p), where the subscripts 

represent the reaction outcome and p is the initial particle concentration, to the 
biased-hopping case. Similar 'duality' relations were found in other formulations 
of A + A — > A or models with microscopically different dynamical rules, 
although for the general case they only apply asymptotically for large times. 

For the 'sticky' two-species reaction, we were able to derive several exact results. 
The concentration decay, for equal initial concentrations of both species, is ~ i -1 / 4 
as in the 'standard' AB case. Furthermore, the exponent value j is not affected 
here by the hopping-rate bias as it is for the 'driven' hard-core AB model. For 
unequal initial concentration of the two species, the asymptotic approach to the 
constant concentration difference is t _3//2 . This clearly contrasts with the expo- 
nential dependence of the 'standard' AB model. Thus, for the multiple-occupancy, 
sticky-particle 'bosonic' model the effects of the hopping bias are not fundamentally 
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important. In fact, in the large-time asymptotic regime the changes due to bias can 
be fully absorbed in the diffusion constant value. This behavior is quite different 
from that of the single-occupancy, hard-core 'fermionic' models. 

To summarize, we have reviewed some common diffusion-limited reactions with 
emphasis on exact results and our recent exact solution of the 'sticky' two-species 
model. Many open problems challenge researchers in this field. The few exact results 
available yield insight on the physical problem and provide a guide for approximate 
analytical and numerical work. We have emphasized the differences between the 
three two-species reactions considered which actually belong to different universality 
classes. 
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